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pO. ABSTRACT 

A computer  simulation  study  of  the  excitation  of  lower  hybrid  waves  by  external 
sources  and  the  associated  plasma  heating  is  presented.  A plasma  slab  with  a 
nonuniform  density  profile  is  modeled  with  a 2 1/2  dimensional  electrostatic 
particle  code.  Both  a finite-length  electrostatic  antenna  and  a single  wave 
exciter  are  considered.  The  frequency  is  chosen  so  that  a wave  conversion  layer 
exists  in  the  plasma.  It  is  found  that  a narrow  sheath  is  formed  at  the  plasma 
edge  and  the  dominant  effect  is  surface  electron  heating  within  that  region. 
Ponderomotive  force  effects  as  well  as  the  nonlinear  generation  of  harmonics  and 

subharmonics  also  occur.  A low  amnlitude  lower  hybrid  wave  is  observed  to  be 
pyritpH  hy  umaom  and  propagates  to  the  wave  conversion  lavpr  where  it  i 
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Abstract 

V 

A computer  simulation  study  of  the  excitation  of  lower  hybrid 
waves  by  external  sources  and  the  associated  plasma  heating  is  presented. 

A plasma  slab  with  a nonuniform  density  profile  is  modeled  with  a 2 1/2 
dimensional  electrostatic  particle  code.  Both  a finite-length  electro- 
static antenna  and  a single  wave  exciter  are  considered.  The  frequency  is 
chosen  so  that  a wave  conversion  layer  exists  in  the  plasma.  It  is  found 
that  a narrow  sheath  is  formed  at  the  plasma  edge  and  the  dominant  effect 
^s  surface  electron  heating  within  that  region.  Ponderomotive  force  effects* 
as  well  as  the  nonlinear  generation  of  harmonics  and  subharmonics  also  occur. 
A low  amplitude  lower  hybrid  wave  is  observed  to  be  excited  by  the  sources 
and  propagates  to  the  wave  conversion  layer  where  it  is  absorbed. 
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I • Introduction 


At  the  present  time  thert  exists  plans  throughout  the  world  for 
high  power  plasma  heating  experiments  using  lower  hybrid  waves.  This 
heating  scheme  has  been  extensively  studied  for  more  than  a decade,  and 
the  present  consensus  is  that  this  is  a difficult  problem,  both  from  the 
experimental  and  theoretical  standpoints.  The  problem  involves  strong 
surface  nonlinearities  due  to  the  large  electric  fields  expected,  and  the 
complications  of  device  geometry  and  plasma  inhomogeneities  make  the 
study  of  wave  propagation  rather  complex  even  in  the  linear  regime. 

Computer  simulation  can  be  a useful  tool  in  understanding  some 
physics  details  that  are  part  of  the  overall  scheme,  especially  those 
aspects  which  are  difficult  to  treat  analytically.  However,  to  do  this 
economically  with  present  day  computers  one  must  idealize  the  geometry 
and  reduce  the  scale  lengths  to  manageable  proportions.  Simulation  thus 
provides  a restricted  but  detailed  view  of  the  problem.  In  the  work  pre- 
sented in  this  paper,  the  computer  model  calculates  the  self-consistent 
particle  trajectories  of  electrons  and  ions  in  a magnetized  slab  geometry 
and  provides  a complete  picture  of  the  dynamics  over  a short  time  scale. 

In  a previous  paper, * we  examined  the  excitation  of  lower  hybrid  waves 
and  their  propagation  in  a homogeneous  slab  and  focused  our  attention  on 
the  coupling  to  bounded  plasma  resources.  For  the  lower  hybrid  heating 
scheme  however,  the  plasma  density  nonuniformity  plays  an  essential  role. 

In  general,  the  density  nonuniformity  gives  rise  to  much  richer  phenomena 
than  is  allowed  in  the  uniform  case  previously  studied.*  In  particular,  one  is 
interested  in  the  case  in  which  the  plasma  contains  a resonance  layer  foT  the 
RF  frequency  used.  The  study  presented  here  is  therefore  a step  closer  to  the 
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experimentally  relevant  situation.  One  other  particle  simulation  on 

2 

this  topic  has  been  reported.  However,  some  of  their  heating  results 
are  different  from  those  in  the  present  study.  For  the  plasma  parameters 
used  in  that  work  = 3,  and  thus  direct  ion  cyclotron  damping  was 

important . 

The  paper  is  organized  as  follows:  Section  II  contains  a discussion 
of  the  computer  model  and  of  the  choice  of  parameters  used  in  the  simulation. 
In  Sec.  II,  we  discuss  the  wave  excitation  and  plasma  heating  with  a finite 
length  source.  The  main  effect  observed  is  that  most  of  the  energy  is  de- 
posited at  the  plasma  surface  and  that  a broad  spectrum  of  harmonics  is 
nonlinearly  generated.  A low  amplitude  lower  hybrid  wave  is  observed  to 
propagate  to  the  resonance  layer  and  is  absorbed.  In  Sec.  IV,  the  case  of 
a single  k exciter  is  discussed.  Finally,  a summary  is  contained  in  Sec.  V. 
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II.  Computer  Model 


An  electrostatic  particle  simulation  model,  described  in  detail 
3 

elsewhere,  has  been  developed  for  this  study  and  pertinent  details  are 
briefly  summarized  here.  The  model  is  bounded  and  is  2Jj  dimensional,  i.e., 
the  field  quantities  and  particle  positions  are  two  dimensional,  while  all 
three  velocity  components  are  retained  in  the  particle  orbits.  The  fields 
are  periodic  in  the  direction  of  the  external  magnetic  field  (z)  and  bounded 
in  the  perpendicular  direction  (x) . Finite-sized  (Gaussian)  particles  are 
used  throughout. 

For  the  simulation  results  presented  here,  vacuum  boundary  conditions 

A 

(i.e.,  no  charges  exist  outside  the  slab  boundary)  are  imposed  in  the  x 
direction  at  x = 0.  By  imposing  symmetry  about  the  plane  of  maximum  density, 
one  needs  to  consider  only  half  the  physical  length  2L^  in  the  x direction, 
reducing  the  computational  effort.  Such  a boundary  condition,  obtained  by 
requiring  the  normal  electric  field  to  vanish  at  x = Lx,  is  "transparent", 
i.e.,  one  avoids  boundary  effects  such  as  sheaths  and  surface  waves.  The 
physical  consequences  of  this  boundary  condition  are  readily  understandable 
in  terms  of  equivalent  image  charges  on  the  other  side  of  the  plane  of 
symmetry.  External  sources,  described  later,  are  specified  at  x = 0,  and 
due  to  symmetry,  there  is  an  equivalent  image  souce  at  x = 2Lx. 

The  choice  of  orientation  of  the  magnetic  field  ]£  , which  is  constant 
in  magnitude,  means  that  while  plasma  drifts  such  as  the  t x f drift  and 
diamagnetic  drifts  do  occur  in  the  model,  their  direction  is  in  the  third 
dimension  (y)  and  therefore  they  cannot  couple  with  any  electric  fields.  This 
configuration  is  useful  for  a first  study,  since  it  allows  one  to  isolate  the 
propagation  of  lower  hybrid  waves  and  plasma  heating  without  the  complications 
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*> 

due  to  parametric  effects  driven  by  E x B drifts  and  scattering  by 
drift  waves,  which  are  not  permitted  to  occur  in  this  geometry. 

In  a future  study,  these  effects  will  be  considered  by  tilting  the 
direction  of  the  magnetic  field  in  the  yz  plane. 

A 

The  initial  plasma  density  is  taken  to  be  uniform  in  the  z 
direction  and  inhomogeneous  in  x,  so  that  ^nABo.  The  density  has  its 
minimum  value  at  the  vacuum  boundary  (x  = 0)  and  increases  monotonically 
to  a maximum  value  at  the  symmetry  plane  x=lx  (with  nmax/nrain  = 10), 
according  to  the  general  form: 


n(x)  = n (B  + A exp  [-  (x  - Lx)2/2a2]), 


where  n is  the  average  plasma  density,  defined  as  the  ratio  of  the  total 

number  of  particles  in  the  simulation  by  the  plasma  area.  The  initial 

temperature  profiles  are  uniform  in  space  with  Tg/Tj  = 10.  The  initial 

electron  thermal  velocity  was  chosen  to  be  v^^  = 0 . 7 <5  w , where  5 is 

the  distance  between  grid  points  used  in  the  simulation,  and  w is  the 

— — 2 % 

average  plasma  frequency,  w = (4^ne  /n>e)  . The  dimensionless  average 
Debye  length  which  will  appear  in  later  expressions  then  has  the  value 


Each  finite-size  particle  has  a distribution  of  charge  given  by: 
- - 2 „ 


P(r.)  = 


J 2**2 


exp  [-  (r  - r.)  /2a  ],  where  a = .85.  The  ratio  of  the  electron 


cyclotron  to  average  plasma  frequency  is  0 e/Upe  = 1*5,  and  the  ion-electron 
mass  ratio  is  NT/me  = 50.  This  then  gives  an  initial  ion  Larmor  radius 

= v^j/H.  = 1.045.  Particles  are  reflected  elastically  at  both  boundaries, 
with  the  full  dynamics  calculated  correctly  to  second  order  in  the  time 

A ~ - 1 

step  flt  = 0.2  “ . The  reflection  of  particles  implies  there  is  no 

pe 


mechanism  for  net  energy  loss  in  the  system.  The  number  of  grid  points 

used  is  L /<$  = 64  and  L /6  = 128. 
x z 

The  external  source  in  these  simulations  is  determined  by  specifying 
an  external  surface  charge  density  oscillating  with  a specified  frequency 
£i>o  on  the  boundary  of  the  slab.  Because  the  magnitude  of  the  charge  density 
is  fixed  externally,  such  a source  corresponds  to  a constant  current  driver. 

The  phase  of  the  oscillator  is  chosen  so  that  the  field  is  zero  at  t = 0 
(i.e.,  a quasi-adiabatic  turn  on). 

Energy  conservation  was  checked  by  calculating  the  work  done  by  the 
external  source  and  subtracting  this  value  from  the  total  plasma  energy. 

In  all  cases,  energy  is  conserved  to  within  3%. 

JII.  Coupling  by  a Short  Capacitor  Plate  Antenna 

In  the  geometry  outlined  in  the  previous  section,  the  lower  hybrid  waves 
of  interest  have  a finite  but  small  wavenumber  k„  along  $ . Since  the  system 
is  periodic  in  z,  the  wavenumbers  in  that  direction  are  quantized  according 
to  k„  = k^  = 2mir/Lz,  where  m is  an  integer  (mode  number).  In  the  perpendicular 
(x)  direction,  the  wavenumber  kx  is  not  constant  because  the  plasma  density 
varies  with  position  according  to  Eq.(l).  In  the  cold  plasma  limit,  electro- 
static lower  hybrid  waves  can  reach  the  resonance  layer  (k,  -v  «°) , where  the 

frequency  oi  of  the  wave  equals  the  lower  hybrid  frequency  id.,,,  defined  by 

LH 


“pe2/!1e2)' 


(2) 


However,  finite  temperature  effects  introduce  an  additional  mode  of  oscillation 
and  instead  of  a resonance,  one  encounters  a wave  conversion  layer  where  two 
roots  of  the  dispersion  relation  coalesce  and  k^  becomes  complex.  The  lower 


hybrid  heating  scheme  is  based  on  the  idea  that  the  mode  converted  wave 
can  be  efficiently  damped  by  plasma  ions. 

Using  a warm  electron-cold  ion  fluid  model,4  one  can  calculate  the 
position  of  the  wave  conversion  layer  for  the  density  profile  used,  Eq.  (1). 

The  result  of  the  calculation  is  shown  by  the  curves  labeled  (m  = 1,  2) 
in  Fig.  1.  Each  of  these  curves  represents  the  location  of  the  layer  for 
the  mode  number  m as  a function  of  frequency.  One  finds  that  com  is  bounded 
by  iojjj  < iom  < (Wp^2  + fl^2)  2,  and  these  upper  and  lower  bounds,  along  with 
the  density  profile  n(x)  are  also  shown  in  Fig.  1.  The  frequency  of  the 
source  u)o  is  chosen  for  practical  reasons  to  be  uo/tope  = 0.14,  which  gives 
a wave  conversion  layer  for  mode  m = 1 about  2/3  of  the  way  inside  the  plasma, 
as  indicated  by  the  dashed  line  in  Fig.  1. 

Two  capacitor  plates  of  finite  length,  oscillating  out  of  phase,  are 
chosen  as  the  coupler  structure.  The  vacuum  equipotential  lines  generated 
by  this  source  are  shown  in  Fig.  2(a).  The  wavenumber  spectrum  of  this  source, 
illustrated  in  Fig.  2(b),  is  a broad-band  spectrum  peaked  at  mode  m = 3 (with 
parallel  phase  velocity  w0Am  = 1-4  v ^ ) . The  magnitude  of  the  surface  charge 
density  a on  each  plate  is  fixed  to  be  eo/kTe  = 0.5761.  This  source  creates 
an  electric  field  at  the  boundary  whose  maximum  amplitude  is  (in  vacuum) 


e 

When  integrated  over  the  plasma  volume,  the  total  field  energy  of  the  source 
(in  vacuum)  is  close  to  7%  of  the  total  electron  kinetic  energy  at  t = 0. 

The  system  is  followed  for  about  11  oscillation  periods  of  the  source. 
The  total  electric  field  energy  (w  ^) , shown  in  Fig.  3,  shows  no  growth  in 


nil 

87TfllcT 


time,  as  when  one  excites  a bounded  plasma  resonance.  The  energies  in  the 
simulation  are  normalized  to  the  initial  total  electron  thermal  energy,  i.e.. 
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kT 

2 p e 


(4) 


for  the  electric  field  energy,  where  = 17280  is  the  number  of  electrons 
in  the  simulation.  The  kinetic  energy  of  species  a is  normalized  in  an 
analogous  manner: 

1 l 


Wa  ~ 3N  kT 
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(5) 


The  electron  kinetic  energy  (wg) , shown  in  Fig.  4,  increases  rapidly  at 
first,  and  then  gTows  at  a slower  rate,  while  the  ion  kinetic  energy  (w^)  shows 
a secular  increase.  At  the  end  of  the  simulation,  the  total  plasma  energy 
increased  by  35%,  with  55%  of  this  energy  gain  going  to  the  electrons  and  45% 
to  the  ions.  The  details  of  this  energy  partition  are  discussed  in  the  next 
section. 

Because  lower  hybrid  waves  were  only  weakly  excited,  two-dimensional 
interferograms  were  made  to  observe  the  space-time  evolution  of  the  plasma 
response  clearly.  These  are  essentially  correlations  between  the  external 
source  and  the  plasma  potential  4>(x,y)  defined  by: 

1 fT 

'Kx.y.f)  - y J0  sin  (wQt)  $ (x.y.t  + f)dt.  (6) 


The  length  of  integration  T must  be  long  enough  to  resolve  the  frequency 
(*iQ  , but  not  so  long  that  slow  changes  in  the  zero  order  plasma  properties 
are  included.  Such  an  interferogram  is  shown  in  Fig.  5,  for  r = 0.  The 
most  striking  feature  in  this  diagram  is  the  large  potential  drop  at  the 


edge,  i.e.,  the  plasma  forms  a sheath  which  shields  out  the  external 
fields.  A much  less  evident  feature  (due  to  the  scale  normalization) 
is  the  low  amplitude  lower  hybrid  wave  which  is  being  excited.  It  is 
more  easily  seen  in  the  equipotential  plot  inserted  in  Fig.  5.  The  well- 
known  resonance  cones  predicted  from  cold  plasma  theory  were  not  observed 
in  this  case  (although  they  have  been  seen  in  our  earlier  studies*  ),  be- 
cause only  a single  mode  penetrates  deep  into  the  plasma.  A Fourier 
analysis  of  the  pattern  shown  in  Fig.  5 in  the  periodic  direction  shows 
(Fig.  6)  that  mode  m = 1 is  the  only  one  which  significantly  penetrates 
into  the  plasma,  because  its  phase  velocity  (uQ/km  =4.1  v ^ ) is  the 
only  one  large  enough  that  it  is  not  subject  to  strong  electron  Landau 
damping.  Although  the  source  had  a broad-band  spectrum  of  modes  (Fig.  2b), 
the  plasma  surface  acts  as  a "filter"  which  transforms  the  spectrum  into 
a single  mode.  This  single  mode  propagates  to  the  wave  conversion  layer, 
which  is  indicated  by  the  arrow  in  Fig.  6. 

From  these  observations  one  concludes  that  there  are  two  distinct 
regions  where  different  plasma  behavior  takes  place.  First  there  is  the 
sheath  region  around  the  source,  where  large  electric  fields  are  present 
and  non-linear  phenomena  are  expected  to  be  important.  Secondly,  in  the 
plasma  interior,  there  is  a low  amplitude  lower  hybrid  wave  which  propagates 
toward  a mode  conversion  layer,  where  the  plasma  properties  remain  essen- 
tially linear.  In  the  following  subsections  we  examine  the  physics  of  these 
two  regions  separately. 


A.  Sheath  region 


The  dominant  effect  which  occurs  in  this  region  (0  < x 5 16) 

is  large  electron  surface  heating,  as  is  expected  due  to  the  large  amplitude 

of  the  electric  fields  and  the  shape  of  the  spectrum.  The  time  evolution  of 

the  electron  velocity  distribution  function  in  the  parallel  direction  f (vM) 

in  this  region  first  shows  electrons  quickly  being  pulled  into  the  tails  of 

the  distribution  within  1-j  oscillation  periods  of  the  source.  The  total 

electron  kinetic  energy  (Fig.  4)  shows  a corresponding  rapid  rise.  On  a 

longer  time  scale,  particles  initially  in  the  bulk  of  the  distribution, 

1 dwe  -4 

also  gain  energy,  but  now  at  a linear  rate  given  by  ~ 2.7  * 10  , 

and  the  energy  shown  in  Fig.  4 reflects  this.  By  the  end  of  the  simulation, 
the  plasma  distribution  function  for  electrons  near  the  edge  (Fig.  7)  is 
approximately  a bi-Maxwellian,  where  the  majority  (-60%)  of  the  electrons 
have  an  effective  temperature  which  is  approximately  ten  times  hotter  than 
the  temperature  of  the  background  (40%) . Figure  8 shows  the  spatial  dependence 
(perpendicular  to  $ ) of  the  electron  kinetic  energy  w , normalized  to  the 
local  plasma  density  ng;  clearly,  the  electron  energy  gain  is  localized  to  the 
edge  region. 

Except  for  the  initial  transients,  the  heating  process  here  is  best 

5-7 

described  by  the  term  transit-time  heating.  When  electrons  transist 

through  a localized  region  of  oscillating  electric  field  with  random  phases 
relative  to  the  rf,  some  particles  gain  a velocity  increment  Ay,  while  the 
velocity  of  an  equal  number  is  decreased  by  Av.  However,  the  energy  gain  of 
particles  with  Av  > 0 is  larger  than  the  energy  loss  of  those  with  Av  < 0, 
and  thus  the  energy  of  a distribution  of  particles  will  on  the  average  increase. 
The  theoretical  rate  of  energy  increase  for  a 1-D  system,  where  Av  <<  v 
has  been  shown  to  be  linear  in  time,  which  is  in  agreement  with  the 

simulation  results.  However,  the  parameters  in  the  simulation  were  such 


i 


that  at  the  surface  Av 


max  mu 


5 vt^e  , which  is  not  the  regime  assumed 


by  the  theoretical  expressions. 


Due  to  the  presence  of  large  localized  electric  fields  at  the  surface, 

one  expects  that  ponderomotive  force  effects  should  play  a role.  This  effect 

is  evident  in  Fig.  9,  where  contour  lines  of  constant  density  are  plotted. 

(The  density  is  time-averaged  over  the  interval  100  <w  t < 500) . One  can 

pe 

see  that  near  the  antenna,  the  plasma  has  been  expelled  and  a density  cavity 
formed  (6ne/ng  < 30%) . The  ions  also  increase  their  kinetic  energy  in  the 
region  near  the  antenna,  due  to  the  direct  acceleration  by  the  source  field. 
Ions  which  acquire  a large  perpendicular  velocity  now  have  large  Larmor 
radii,  and  they  migrate  into  the  plasma  interior.  This  surface  ion  heating 
will  be  discussed  in  more  detail  in  a later  section. 

In  summary,  one  attains  a plasma  at  the  edge  whose  electron  temperature 
is  ten  times  hotter  and  in  which  a density  cavity  forms  around  the  antenna  due 
to  the  ponderomotive  force.  These  changes  at  the  edge  then  modify  the  manner 
in  which  the  ’deal  capacitor  plate  antenna  couples  the  external  RF  energy  to 
the  plasma  interior. 

B.  Spectral  Broadening 

Another  important  feature  which  has  its  origin  in  the  sheath  region 
is  the  generation  of  a broad  frequency  spectrum  which  includes  harmonics  and 
subharmonics  of  the  external  excitation  frequency  In  Fig.  10,  we  show  a 

power  spectrum  for  mode  m = 1 at  various  points  in  space,  defined  by 


P_(x,w)  = 


4it  n kT 
o e 


\K  (x,m)|2- 


(7) 


is  the  dominant  feature  in  the 


One  finds  that  the  driving  frequency 

spectrum  for  those  values  of  x such  that  w < id...  (x),  i.e.,  below  its 

cut-off  point  at  the  wave  conversion  layer  at  x = 40.  In  addition,  one 

finds  ion  cyclotron  modulation  which  is  evident  in  the  generation  of 

upper  and  lower  side-bands  w » to  ± fi. . Since  the  lower  side-band  is  cut 

o 1 

off  at  a location  given  by  to  - fi.  * toIU  (x)  < to  (i.e.,  at  x = 24),  it 

does  not  penetrate  deep  into  the  interior;  however,  the  upper  sideband 

penetrates  deeper  than  the  fundamental  wave  launched  by  the  antenna 

(e.g.,  at  x = 48).  This  ion  cyclotron  modulation  is  a phenomena  which  we 

observed  in  our  earlier  study  whenever  finite-size  antennae  were  used. 

Another  frequency  seen  near  the  sheath  region  (x  < 12)  in 

Fig.  10  is  the  half-harmonic  w /2.  It  turns  out  that  the  local  lower 

o 

hybrid  frequency  near  the  antenna  is  approximately  equal  to  wq/2.  A 
plasma  oscillation  at  u)q/2  is  found  to  be  excited  and  driven  to  large 
amplitude  near  the  density  depression  created  by  the  ponderomotive  force 
(Fig.  9).  This  excitation  process  is  an  example  of  the  classic  parametric 
excitation  described  by  Mathieu's  equation: 


+ n 2[1  + e cos  u tlx  = 0 

dt2  0 ° 


(8) 


when  a normal  mode  frequency  is  modulated  by  an  external  source.  3 

This  model  equation  predicts  that  when  u>o  = 2Uq  odd  numbered  harmonics  of 

u)q/2  would  also  be  excited,  e.g.,  3wo/2  arises  due  to  the  non-linear  mixing 


of  a)0/2  with  the  fundamental  u)Q.  In  Fig.  10,  one  indeed  sees  peaks  at 
3u)q/2,  5<do/2  and  7^/2.  Finally,  the  first  harmonic  at  2Wq  is  also  seen 
in  Fig.  10,  corresponding  to  the  non-linear  mixing  of  the  fundamental 
with  itself.  One  can  estimate  from  Fig.  10  that  the  power  in  the  harmonics 


is  an  order  of  magnitude  smaller  than  the  power  at  u)0- 

Harmonics  at  other  mode  numbers  (m  = 0,  2)  were  also  observed.  For 
mode  m = 2,  we  find  that  a peak  at  w = Wo  exists  only  near  the  sheath  region 
because  <jj0/k,,  = 2.0  for  this  mode,  as  expected  from  Fig.  6.  However, 

the  first  harmonic  w = 2ojq  is  found  in  the  plasma  interior  because  it  has 
a phase  velocity  (2^/k,,  = 4.1  vt^  ) in  excess  of  the  electron  thermal  velocity 


and  is  not  subject  to  strong  Landau  damping.  The  frequencies  w0/2,  3U  / 2 and 
5o>o/2  were  also  observed  for  this  mode. 

C.  Wave  Propagation  in  the  Plasma  Interior 

Wave  propatation  characteristics  can  be  extracted  by  introducing 
a lag  time  x in  the  interferogram  defined  by  Eq.  6.  In  Fig.  11,  we  examine 
the  propagation  of  the  dominant  mode  m = 1,  by  showing  the  spatial  dependence 
of  the  potential  at  succeeding  lag  times.  One  can  see  that  the  phase  fronts 
are  propagating  toward  the  antenna.  Since  lower  hybrid  waves  have  oppositely 


directed  group  and  phase  velocities,  this  confirms  that  the  direction  of 
energy  flow  is  toward  the  interior.  The  wave  propagates  toward  the  wave 
conversion  layer  with  diminishing  amplitude  until  it  is  totally  absorbed. 

To  examine  more  carefully  whether  any  wave  conversion  to  short  wavelength 
ion  waves  occurs  as  the  wave  approaches  the  center  of  the  plasma,  we  show  in 
Fig.  12,  the  interferograms  of  the  ion  density  (x,x),  defined  in  analogy 
to  Eq.  6.  First,  one  sees  that  the  ions  contribute  to  the  lower  hybrid  wave 
only  when  the  density  is  sufficiently  large  (x  - 30,  o)q  < wLH) • When  the 

| 
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density  is  low,  the  lower  hybrid  wave  is  predominantly  an  electron  wave. 

In  the  high  density  regions,  where  the  ions  are  contributing  to  the  lower 
hybrid  wave,  one  finds  that  the  phase  fronts  are  moving  in  the  -x  direction, 
with  the  same  wavelength  as  in  Fig.  11.  Since  no  appreciable  standing  wave 
pattern  is  observed,  one  can  conclude  (to  within  10%)  that  neither  reflection, 
nor  transmission  from  the  image  source  at  x = 2Lx  (due  to  the  use  of  symmetric 
boundary  conditions)  has  occurred.  The  conclusion  is  that  the  wave  is  totally 
absorbed  with  no  measurable  wave  conversion. 

In  spite  of  the  fact  that  the  wave  is  totally  absorbed,  it  is  difficult 
to  show  any  corresponding  increase  in  particle  kinetic  energy,  since  the  wave 
amplitude  is  so  low  (only  a few  times  the  thermal  noise  level)  and  other  effects, 
such  as  heat  exchange  due  to  collisions  between  the  initially  hot  electrons  and 
the  cold  ions,  are  more  important. 

IV.  Lower  Hybrid  Waves  Driven  by  Single  k Exciter 

In  order  to  excite  the  lower  hybrid  wave  more  efficiently,  we  chose  an 
external  source  where  the  energy  would  be  predominantly  in  the  mode  m = 1. 

The  geometry  consisted  of  two  capacitor  plates  each  of  length  Lz/2.  The 
amplitude  was  chosen  (ea/xTe  = .186  *)  so  that  the  total  energy  of  the  source 
(in  vacuum)  was  the  same  as  the  previous  case,  but  concentrated  primarily  in  a 
single  mode.  All  the  other  simulation  parameters  remained  the  same. 

The  two  dimensional  interferogram  (Fig.  13)  for  this  case  shows  that 
although  the  structure  of  the  sheath  has  changed,  the  low  amplitude  lower 
hybrid  wave  excited  in  the  plasma  interior  looks  much  the  same  as  before,  with 
about  a factor  of  two  increase  in  amplitude.  The  major  difference  from  the 
previous  case  was  in  the  properties  of  the  sheath  region. 


The  kinetic  energy  of  each  species  as  a function  of  time  is  shown 
in  Fig.  14.  One  can  see  that  while  the  electron  kinetic  energy  gain  was 
more  than  double  the  previous  case,  the  ion  energy  gain  was  about  the  same. 

As  mentioned  previously,  the  ions  gain  energy  by  being  accelerated  directly 
by  the  source  fields.  This  effect  is  illustrated  in  the  phase-space  dia- 
gram in  Fig.  15,  where  we  plot  the  co-ordinates  (x,  v^)  for  all  the  ions. 

One  finds  a small  population  of  ions  being  accelerated  from  the  source  at 
x = 0 to  velocities  as  large  as  fifteen  times  the  ion  thermal  velocity,  and 
penetrating  to  the  plasma  interior.  (In  a real  experiment,  of  course,  many 
of  the  ions  would  be  lost  in  the  system,  whereas  in  this  computer  experiment 
the  ions  are  reflected  if  they  reach  the  boundary).  Figure  16  shows  the 
spatial  distribution  of  the  ion  kinetic  energy  density  in  this  case. 

There  were  fewer  harmonics  generated  with  this  source.  For  mode  m = 1, 
the  ion  cyclotron  modulation  was  absent,  and  the  u>0/2  excitation  was  much 
weaker.  The  first  harmonic  2^  was  also  present.  For  mode  m = 2,  the 
main  effect  was  an  oscillation  with  (2Wo>  2kQ) , since  no  power  was  being 
supplied  by  the  source  at  that  wave  number.  There  is  thus  an  effective 
up-shifting  of  the  wavenumber  spectrum  due  to  the  non-linearities  in  the 
sheath. 

The  wave  propagation  characteristics  of  the  fundamental  in  this  case  were 
the  same  as  those  shown  in  Fig.  11. 

V.  Conclusion 

In  all  RF  heating  experiments  one  must  contend  with  sheath  regions  around 
the  antenna.  Large  amplitude  sheaths  are  especially  difficult  to  analyze  and 
are  poorly  understood.  In  this  simulation,  we  have  demonstrated  that  for  electro- 
static excitation,  such  sheaths  play  an  important  role  in  the  efficiency  of  excita 
tion  and  the  power  absorption.  Further  simulation  can  shed  much  light  on  this 
complex  problem. 


To  summarize  the  results,  with  electrostatic  antennae  of  the  type  used 


here,  most  of  the  energy  lies  within  a narrow  sheath  region  at  the  edge. 

Only  modes  with  u)o/k„  >2.5  v^g  are  excited  in  the  plasma  interior,  and  these 

have  low  amplitudes.  The  result  is  that  the  dominant  effect  is  surface  heating. 

The  electrons  at  the  surface  absorb  55-75%  of  the  total  energy.  Ions  absorb 

most  of  the  remainder  at  the  surface  and  migrate  toward  the  interior.  For 

finite-sized  antennae,  a significant  density  depression  (6n/n  Z 30%)  is  created 

in  front  of  the  source.  This  depression  plays  a role  in  some  of  the  harmonics 

which  are  non-linearly  generated.  About  5%  of  the  energy  goes  into  a wide 

variety  of  harmonics,  e.g.,  wq/2,  3u)q/2,  2u>q.  Ion  cyclotron  modulation  at 

a)  - oj  ± fl.  is  observed  whenever  finite-size  antenna  were  used.  The  low 
o 1 

amplitude  fundamental  wave  at  (dq  is  observed  to  propagate  toward  the 
interior,  where  it  is  absorbed. 


w 


1.  V.  K.  Decyk,  J.  M.  Dawson,  and  G.  J.  Morales,  Phys.  Fluids  22,  507  (1979) 

2.  Hirotada  Abe,  Hiroyuki  Kajitani,  and  Ryohei  Itatani,  Nagoya  University 


Institute  of  Plasma  Physics  Report  IPPJ-296  (1977). 

3.  Viktor  K.  Decyk  and  John  M.  Dawson,  J.  Computational  Phys.,  to  be  published 

4.  Viktor  K.  Decyk,  Ph.D.  dissertation.  University  of  California,  Los  Angeles 

(1977). 

5.  G.  J.  Morales  and  Y.  C.  Lee,  Phys.  Rev.  Lett.  33,  1534  (1974) 

6.  Hideo  Sugai,  Kazuo  Ido,  and  Susumu  Takeda,  J.  Phys.  Soc.  Japan  46,  228  (1979) 

7.  J.  M.  Berger,  W.  A.  Newcomb,  J.  M.  Dawson,  E.  A.  Frieman,  R.  M.  Kulsrud, 

and  A.  Lenard,  Phys.  Fludis  1_,  301  (1958). 

8.  Mechanics  by  L.  D.  Landau  and  E.  M.  Lifshitz  [Addison-Wesley,  Reading, 

Mass.,  1969],  p.80. 

and  C.  S.  Liu,  Advances  in  Plasma  Physics,  vol  6, 

New  York,  1976],  p.7. 


9.  Kyoji  Nishikawa 
[John  Wiley, 


j 


i 


FIGURE  CAPTIONS 


Fig.  1 


Lower  Hybrid  wave  conversion  frequencies  u>  (x)  for  mode  numbers 

111 

m = 1,  2,  the  frequencies  a)T„(x)  and  /ai  .^(x)  + , for  the 

Ln  pi  l 

initial  density  profile  n (x) . X is  the  co-ordinate  perpendicular 
to  the  magnetic  field.  The  driving  frequency  (oo  is  chosen  to 
give  a wave  conversion  layer  for  m = 1 about  2/3  of  the  way  inside 
the  plasma. 


Fig.  2 (a)  Equipotential  lines  for  finite-size  exciter  in  vacuum. 

(b)  Wave  number  spectrum  for  finite-size  exciter  in  vacuum.  The 
small  secondary  peaks  for  large  wavenumber  are  due  to  the  sharp 
cutoff  of  the  external  charge  distribution  on  the  plates. 

Fig.  3 Time  evolution  of  the  total  electric  field  energy  for  the  finite- 
size  exciter  case. 


Fig.  4 Time  evolution  of  the  electron  (we)  and  ion  (w^)  kinetic  energies 
for  the  finite-size  exciter  case.  (Initial  temperature  ratio  was 
T /T.  = 10.) 

Fig.  5 Spatial  dependence  of  the  interferograms  of  the  plasma  potential 
(at  x = 0)  with  equipotential  plot  inserted,  showing  the  large 
sheath  region  around  the  antenna  and  the  low  amplitude  lower  hybrid 
waves,  for  finite-size  exciter  case. 

Fig.  6 Interferograms  of  plasma  potential  (at  t = 0)  for  Fourier  modes 
m = 1,2,  shown  as  a function  of  the  perpendicular  co-ordinate  X. 

(Xn  * .7;  solid  line  shows  Im  4>  and  dashed  line  shows  Re  <(>  ) . 


Electron  velocity  distribution  parallel  to  the  magnetic  field  for 
electrons  near  the  edge  region  (0  < x < 16).  (Dashed  line  shows 
initial  distribution  and  solid  line  shows  distribution  at 


Fig.  7 


to  t = 480)  . 
pe 

Fig.  8 Spatial  dependence  perpendicular  to  Bq  of  the  electron  kinetic 
energy  density  at  various  times. 

Fig.  9 Contour  lines  of  constant  values  of  the  time  averaged  electron  density 
ng,  illustrating  the  density  cavity  formed  around  the  antenna. 

Fig.  10  Power  spectrum  for  Fourier  mode  m = 1 as  a function  of  frequency 
at  various  spatial  locations.  <oq  is  the  externally  excited  mode. 

Fig.  11  Interferograms  of  potential  due  to  Fourier  mode  m = 1 as  a function 
of  the  perpendicular  co-ordinate  at  various  lag  times  x shows  direc- 
tion of  wave-phase  propagation  corresponds  to  a backward  wave. 

Fig.  12  Interferograms  of  ion  density  for  Fourier  mode  m = 1 as  a function 
of  the  perpendicular  co-ordinate  at  various  lag  times  x. 

Fig.  13  Spatial  dependence  of  the  interferograms  of  the  potential  (at  x = 0) 
with  equipotential  plot  inserted,  for  single  k exciter. 

Fig.  14  Time  evolution  of  the  electron  (w  ) and  ion  (w.)  kinetic  energies 

6 1 

for  the  single  k exciter  case.  (Initial  temperature  ratio  was 
Te/T.  = 10). 

Fig.  15  Ion  phase  space,  showing  the  co-ordinates  (x,  v ) of  all  the  ions  at 

time  w t = 100. 
pe 

Fig.  16  Spatial  dependence  perpendicular  to  Bq  of  the  ion  kinetic  energy 
density  at  various  times. 
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